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ABSTRACT 

The study of the evolution of a suprathermal electron beam traveling through 
a background plasma is relevant for the physics of solar flares and their associ- 
ated type III solar radio bursts. As they evolve guided by the coronal magnetic 
field-lines, these beams generate Langmuir turbulence. The beam-generated tur- 
bulence is in turn responsible for the emission of radio photons at the second 
harmonic of the local plasma frequency, which are observed during type III solar 
radio bursts. To generate the radio emission, the beam-aligned Langmuir waves 
must coalesce, and therefore a process capable of re-directioning the turbulence 
in an effective fashion is required. Different theoretical models identify the elec- 
trostatic (ES) decay process Li ^ L2 + S (L: Langmuir wave; S: Ion-acoustic 
wave) as the re-directioning mechanism for the L waves. Two different regimes 
have been proposed to play a key role: the back-scattering and the diffusive 
(small angle) scattering. This paper is a comparative analysis of the decay rate 
of the ES decay for each regime, and of the different observable characteristics 
that are expected for the resulting ion-acoustic waves. 

Subject headings: Sun: corona — Turbulence — Sun: radio radiation 



1. Introduction 

During solar flares, large amounts of energy are released and transformed in coronal 
heating and particle acceleration, where regions of magnetic reconnection are believed to 
be the acceleration sites for suprathermal electron beams. Once accelerated, these beams 
travel through the coronal plasma guided by the coronal magnetic field-lines, and generating 
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a variety of observable emissions. One long standing discrepancy in the modeling of this 
phenomena is that a beam with a given energy flux seems not to be able to simultaneously 
reproduce the observed emissivities in HXR (due to non-thermal bremsstrahlung at the chro- 
mosphere) and radio emission (due to beam-generated Langmuir turbulence in the corona). 
Once the beam energy flux is set to reproduce the HXR levels, the derived radio emissivities 
tend to be much higher than observed levels (Emslie & Smith 1984; Hamilton & Petrosian 
1987). 

In this context, we have developed a model for the evolution of electron beams and 
the generation of Langmuir turbulence, and we have also computed the emission of radio 
waves due to the coalescence of the beam-excited Langmuir waves. Our models describe the 
evolution of the beam and the produced Langmuir turbulence, consistently considering the 
effect of both collisions and quasi-linear wave-particle interaction. The level of turbulence 
derived from our models is up to two orders of magnitude lower than previous attempts 
(Vasquez & Gomez 1997). The production of photons at the second harmonic of the plasma 
frequency (radio waves), is the result of the coalescense ot two Langmuir waves. In our 
model we assume that the beam-generated Langmuir waves become isotropic in an effective 
way, due to electrostatic decay Li + L2 — > T(2c<jpc). We found that an adequate treatment 
of the second harmonic photons generation (relaxing the head-on approximation) yields to 
further reductions of the radio emission (Vasquez et al. 2002). Although these results help 
to reduce the gap between the predicted HXR and radio emissivities, we find that further 
reductions are required. 

The observed radio emission requires the coalescense of the beam-generated Langmuir 
waves. Therefore, a process capable of re-directioning the turbulence in an effective fashion 
is required. Different models in the hterature resort to the electrostatic (ES) decay Li — > 
L2 + S {L: Langmuir wave; S: Ion-acoustic wave) as the re-directioning mechanism for 
the L waves. Two different regimes have been proposed to play a key role. One of them 
is the so called back-scattering limit of the ES decay (Edney & Robinson 2001; Cairns 
1987). In this limit, the primary Langmuir wave decays into another one that propagates 
almost in the opposite direction. The other asymptotic regime is the small-angle limit of 
the ES decay, which has also been pointed out by some authors as potentially relevant 
in this context (Melrose 1982; Tsytovich 1970). In this hmit, the propagation directions 
of both Langmuir waves (Li and L2) form a small angle, and the ES decay acts on the 
beam-generated Langmuir waves as a difussion mechanism through k-space. In our models 
described above, the isotropization of Langmuir waves has been an assumption rather than 
a result obtained from our beam-turbulence quasi-linear equations. If the timescale of the 
ES decay is short enough (as compared to the timescale for the Langmuir waves to escape 
their generation region), it is reasonable to expect that the small-angle limit will render the 
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beam-generated Langmuir turbulence isotropic. In this work we revise this approximation 
in detail. We compare the rate of ion-acoustic wave generation in both limits, for a given 
set of beam-generated Langmuir spectra. We also analize the resulting frequencies of the 
ion-acoustic waves produced in each limit, and compare against reported in-situ observations. 

In-situ observations of type III solar radio bursts have shown clear evidence in support 
of the ocurrence of the ES decay simultaneously with Langmuir waves. For example. Cairns 
& Robinson (1995) have analyzed ISEE 3 data that shows the coexistence of high and low 
frequency electrostatic waves, identified as Langmuir and ion-acoustic waves respectively, 
during a type 111 radio event. They analyze the observed low frequency waves, and find 
that their frequencies are consistent with those predicted by assuming the existence of the 
ES decay acting in the back-scattering limit. Similar observational works by Thejappa et 
al. (2003) and Thejappa & MacDowall (1998), analyzing URPE/ Ulysses data, also support 
the occurence of the ES decay in association to impulsive Langmuir waves excited during 
type III radio events. In the present paper wc make a similar analysis to the one by Cairns 
& Robinson (1995), and find out that the URPE data analyzed by Thejappa & MacDowall 
(1998) and Thejappa et al. (2003) is consistent with the occurrence of the ES decay acting 
in the diffusive limit. 

This paper is organized as follows. In section 2 we revisit the ES decay and the general 
expression of its rate of occurrence. In section 3 we apply the general results to the specific 
case of a beam-generated Langmuir spectrum. We make a quantitative comparison between 
the diffusive and back-scattering cases. In section 4 we compute the expected frequencies for 
the generated ion-acustic waves for the particular type III solar radio bursts observed and 
studied by Thejappa & MacDowall (1998) and Thejappa et al. (2003). Finally, in section 5 
we list our main conclusions. 



2. Electrostatic Decay Rate 

The electrostatic decay Li ^ L2 + S (where Li2'- Langmuir; S: Ion-acustic), must 
verify the momentum and energy conservation, so that their wave vectors and frequencies 
satisfy 



where the dispersion relationships are. 
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^^L(ki)~c^pe(^l + ^(A:LADe)'j i ns{ks) ^ kgYs (2) 

Prom these relations, we obtain the modulus and direction of the wave vectors of the resulting 
waves L2 and 5" (see Figure 1), with respect to the wave vector of the initial Langmuir wave 
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The spectral density of the plasmons of type a (with a = L, S in this case), is 
defined such that the volumetric density of plasmons is n"^ = J N^. In terms if these 
spectral densities, the rate of production of ion-acoustic (S) waves is (Tsytovich 1970) 
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where the probability of the decay is 
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where rUe, nip are the electron and proton mass respectively, Vre is the electron thermal 
velocity, and the delta functions express the momentum and energy conservation. 



3. Results for beam-generated Langmuir waves 



We use the expresions of the previous section to compute the decay rate of a beam- 
generated Langmuir spectrum. Hereafter we adopt a given shape for a beam generated 
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Langmuir spectrum. The dependence of the spectrum upon the modulus of wavenumber, 
has been obtained by Vasquez & Gomez (1997). Wc also assume this spectrum to by 
axisymmetric and with a gaussian angular spread about the direction of the beam, and 
analyze the initial ES decay rates. 

The beam-generated Langmuir turbulence is initially aligned with the beam propagation 
axis (the local magnetic field). As it proceeds, we expect the ES decay to re-directionate the 
Langmuir waves. We thus divide our analysis in two cases. In a first place we analyze the 
ES decay of a perfectly colimatcd Lanmguir spectrum. Later on, we analize the same rates 
for a non-colimated spectrum. In both stages of our analysis we compute and compare the 
rates for both difussive and back-scattering cases. 



3.1. Colimated Langmuir spectra 

Let us evaluate the initial rate of production of S waves generated by a Langmuir 
spectrum, which is in turn produced by a perfectly colimated beam (i.e. we take N^^^ = 0), 



T(k5) dt ^ J (27r)6 

l2 



where the integral in has been performed using the momentum conservation, and the 
function K is given by 

ne^l]|m,(27r)6 o^pe . . 

In a previous work (Vasquez & Gomez 1997; Vasquez et al. 2002) we developed a model 
to derive the spectrum of beam-generated Langmuir waves. According to this model, the 
wavenumber of the excited Langmuir waves lays in the range 



h = [kmin, A:max] ~ [6, 12] A;o Vn (10) 

where Tg = Te[K]/10^. Approximating to a constant level for the perfectly colimated spec- 
trum, and calling ez to the direction of the beam. 
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= (27r)3^<5(A;i.)5(A;i,), (11) 

where AA;^ is the spectral width, the result for the initial production rate of ion-acoustic 
waves is 
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where the dimensionless function P{ks, //) is given by 
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/X 4 + {ks/koY + 4(A;5/A;o)(l - 2/x^) 2// 

and is zero otherwise. Figure 2 shows the dimensionless decay rate P{ks, //) given by equation 
(13), as a function of the wave number ks and the angle ^sl between kg and k^,. Here we 
have taken Tg = 0.17 as a characteristic number for a type III solar radio bursts analyzed in 
section 4. This graphic confirms that the cases of maximum probability of occurrence are the 
back-scattering (with ^sl --^0, P ~ 1) and the diffusive scattering (with ^sl ~ 75°, P ~ 4). 
The factor of four in the diffusive decay rate results from the modulation factor in 
equation (13). The propagation direction of the resulting ion-acoustic wave is parallel to the 
beam propagation in the case of back-scattering, and almost perpendicular to the beam in 
the diffusive case (with ks <^ kL)- We therefore find that the decay rate becomes maximum 
for two limiting cases: a) small-angle (diffusive) scattering: a — > 0, and b) back-scattering: 
q; — > TT. 

Under the diffusive scattering approximation, estimates for the characteristic isotropiza- 
tion and energy transfer timescales, indicate that the latter is much larger, implying that 
the isotropization occurs in a quasi-elastic fashion Tsytovich (1970) (see also Vasquez et al. 
2002). The quasi-conservation of the total energy is readily seen from the following argu- 
ment. The ES decay conserves the total number of Langmuir plasmons. Now, given the weak 
dependence of their energy on the wave-number (see equation (2)), we can approximate the 
energy of each plasmon as Upe- Therefore, this decay can not produce a significant change 
in the total Langmuir wave energy. 
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3.2. Non-colimated Langmuir spectra 

In the previous section we found that the ES decay of the ID Langmuir spectrum 
generated by a perfectly cohmated beam is initially dominated by small-angle decay processes 
and, in second place, by back-scattering. We then expect that after a transient stage, 
the Langmuir spectrum is not longer perfectly aligned with the beam direction, and that 
a growing back-scattered spectrum also arises. Nonetheless, we expect that the forward 
spectrum will contain much more energy than the backward spectrum, at least during the 
early stages of the evolution. Under this assumption, we consider now an axisymmetric 
spectrum around the beam direction, and model it through functions A^k — N[6, k) that 
peak in the forward direction and monotonically decrease as ^ — > tt. 

Let us start by analyzing the small angle scattering limit, i.e. ks <tC kz- For the diffusive 
decay we have that the spectral width of the ion-acoustic waves produced is much smaller 
than the spectral width of the primary Langmuir waves. More specifically, we have 
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^kL {ki 

where (ks) — ksmax/'^ and (ki) — (A;min + ^max)/2 indicate mean wavenumber values. As the 
decay proceeds, each ES decay process will produce one ion-acoustic wave, hence 



i.e. the spectral density of the (emitting) Langmuir waves becomes much lower than that 
of the (emitted) ion-acoustic waves (even though the Langmuir energy density can be much 
larger than that of the ion-acoustic waves). Thus, we neglect the term N^N^ in the decay 
rate given by equation (6) and approximate 



where K is given by equation (9) (see also equation (8)). Using the energy conservation, 
equation (17) reduces to the double integral. 
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where we have assumed the Langmuir spectrum iV^ to be symmetric about the beam pro- 
pagation axis ez (i.e. independent of the angle 0), and the values k2, 02 are related to ki, 9i 
through 
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and cos (a) is given by equation (3). We note here that in the back-scattering limit, the 
term NliNl2 can not be neglected. Instead, for back-scattering, we expect that Nl2 <^ Nli 
(at early stages of the evolution), and hence we can neglect the two last terms in equation 
(6). Therefore, the rate for back-scattering can be obtained from equation (17) by simply 
neglecting the second term Ng^, and also taking the limits cos (a) — > — 1 and /i ^ +1. 

To compare these rates against the perfectly colimated case (equation (12)), let us 
multiply and divide the RHS of equation (17) by Ki = {27T)^nL/ Ak^, to obtain 

^ ^-Pi{ks,9s) (20) 



T(k5) To' 

where I/tq is the reference rate given by equation (12), 9s is the angle between the ion 
acoustic wave vector and the beam propagation direction ez, and the dimensionless function 
Piiks,9s) is given by 
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We use this expression to evaluate the decay rate for a spectral model that is stronger in 
the forward direction (but not colimated). Here we neglect the dependence in wavenumber, 
and assume the spectral density as constant within the range [/cmin, ^max] given by equation 
(10). More specifically we consider 
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N{k,9) — ^ -j^ _|_ g(6)-6»c)/A6> ' ^min < k < k^ax (22) 

where the normalization constant C is such that the number density of Langmuir plasmons 
is = J (2^^k- As an example, Figure 3 shows this function for 9c = 25°, = 2°. This 
angular distribution peaks at ^ = and monotonically decreases with increasing 6. The 
mean value is reached at 6c- Thus, for smaller values of this parameter the spectrum is more 
colimated. The parameter measures the half-width over which the function varies from 
75% to 25% of its maximum value. For smaller values of this parameter, the spectrum has 
a sharper angular edge. 

Figure 4 shows the resulting Pi for two sharp-edged angular spectra with A^^ = 2°. 
One spectrum is well concentrated along the beam direction with 6c = 5°, and the other 
corresponds to a wider distribution with 9c — 25°. Our expression for Pi (equation (21)) 
is valid in the diffusive regime ks ^ /cl, so we show results for ks — 0, .., (/ci) /lO. Rates 
monotonically increase with ks, until they saturate at about (ki) /5. For larger values of ks 
the rates decrease with increasing ks- 

Higher rates are obtained for a more colimated beam, and the peak value is reached at 
a larger angle ^5. This behaviour can be understood in terms of our analytic decay analysis 
for the perfectly cohmated case. The peak is expected to occur as the result of the decay of 
Langmuir waves located in regions where the angular gradient of the spectrum is maximum, 
i.e. located around 6c- This is because in this region the difference N{ki,6i) — N{k2,62) in 
equation (21) reaches its largest value. On the other hand, for the characteristic numbers of 
type III solar radio bursts being used in this analysis, our results (sec Figure 2) show that 
the ion-acoustic waves are preferentially directed toward angles around the value 63 ~ 75° 
(with respect to the direction of the decaying Langmuir wave). We can therefore expect 
then that Langmuir waves with an angle dc respect to the ez axis will preferentially produce 
ion-acoustic waves with an angle ds ~ 75° — 6c- This rough estimate gives peak angles of 
70°and 50° for the two values of 6c considered here, which is in agreement with the detailed 
numerical results shown in Figure 4. 

Also, in the perfectly colimated limit {6c — > 0, A^ — > 0), from comparison with equation 
(13), we expect that peak values of Pi will be close to the peak value P ~ 4 corresponding to 
the perfectly colimated case (sec Figure 2). The numerical results of Figure 4 are consistent 
with this. As a consistency check. Figure 5 shows the same results for a highly colimated 
distribution with 6c = 0.5°, A^ = 0.5°, with peak rates of Pi = 3.5, at 6s ^ 74°. 

For a fixed ks — {ki) /lO, Figure 6 shows the effect of reducing the Langmuir wave 
distribution angular gradient. For the distribution with the edge centered on the angle 
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dc — 5°, left panel shows the resulting rates for = 2°, 10°. The right panel shows the 
corresponding results for the distribution with edge on the angle O^. = 25°. The effect of a 
smoother distribution is to extend the possible angles of the produced ion-acoustic waves 
over a much wider range. Also, as the angular gradient of the Langmuir spectrum is reduced, 
the difference Ni — N2 in equation (21) decreases, and hence Pi {9s) yields lower values. 

Let us also compare the different diffusive decay rates against the corresponding back- 
scattering decay rates. The backscattering of a Langmuir plasmon of wavenumber pro- 
duces an ion-acoustic plasmon of wavenumber ks = 2A;^ — 2kQ, as follows from taking the 
limit /X — > -|-1 in equation (4). Figure 7 shows the resulting back-scattering Pi for the same 
two Langmuir spectra analyzed in Figure 4. Thin curves correspond to ks = 2kL,min — 2/co 
and thick curves to ks — ^k^^max ~ 2/co. In this case, being the result of back-scattering, the 
angular distribution of the decay rate closely follows that of the decaying Langmuir waves 
N{ki, 9i) in equation (21) (compare for example the right panel of Figure 7 with Figure 3). 

A comparison between the back-scattering results of Figure 7 against the diffusive decay 
results Figure 4 shows that the diffusive rates are always larger, by factors of order 2 to 5. 
Figure 8 shows the same results for a highly colimated distribution with ^c = 0.5°, A^ = 0.5°. 
The resulting rates are of order Pi ~ 1.5, or about a factor of 2 lower than the corresponding 
diffusive decay values of Figure 5. These results arc consistent with the perfectly colimated 
beam results of Figure 2. Therefore we find that, in all cases, diffusive decay rates are 
systematically larger (though comparable) than the corresponding back-scattering rates. 

4. Frequencies of the ion-acoustic Waves 

Several space-based missions, located at distances of order 1 AU from the Sun, are 

designed to perform in-situ measurements of the extended solar wind. These missions are 
able then to perform in-situ measurements of type 111 solar radio bursts, that sometimes 
extend far away from the Sun and reach the Earth. One of these instruments is the Unified 
Radio and Plasma Wave Experiment (URAP) aboard the Ulysses mission. For the purpose 
of interpreting type 111 burst data collected by this experiment, let us suppose that an 
ion-acoustic wave travels through the region of the experiment. If Vgw is the local solar 
wind velocity, the Doppler-shifted frequency of the ion-acoustic wave, as measured by the 
experiment, is given by 



fs = ^{ksVs + ks-\,^)^ 
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where ks has been ehminated from equation (4), ^sl and ^sr are the angles between the 
propagation direction of the ion-acoustic wave and those of the primary (beam ahgned) 
Langmuir waves, and the (radial) solar wind, respectively. This formula is analogous to 

the one used by Cairns & Robinson (1995) to analyze a type III solar radio burst. They 
assume the ES decay products only in the back-scattering limit, and hence they consider 
the particular case cos (^sl) = +1- Also in that limit, the produced ion-acoustic waves 
propagate aligned with the beam direction, so that 9sr = Olv, i-e. the angle between the 
primary Langmuir waves and the solar wind direction. On the other hand, in the diffusive 
hmit, ^SL > so that the angle ^sr takes the range of values |^sl — ^Lr| < ^Sr < ^sl + ^Lr- 

To quantify the predicted frequencies for the ion-acoustic waves in each of these two 
limiting cases, we refer to a specific observational case, analyzed by Thejappa & MacDowall 
(1998). They show the data of a type III solar radio event recorded by URAP/ Ulysses on 
March 14 1995. Their detailed analysis shows radio emissions at both the plasma and 
second harmonic frequencies. At the same time, impulsive electrostatic fluctuations were 
recorded at the local plasma frequency, which turns out to be /pe ~ 23 kHz, identifled 
as Langmuir wave bursts. The WFA (Wave Form Analyzer) instrument detected electric 
fluctuations, highly correlated in time with the observed Langmuir impulsive peaks, in the 
range 450 Hz. The instrument also detected magnetic fleld fluctuations in the range 
— > 50 Hz. The combination of these facts indicates that the electric fleld fluctuations in the 
range 50 — >■ 450 Hz are of a pure electrostatic nature, such as ion-acoustic waves. If this is the 
case, its strong temporal correlation with the Langmuir bursts supports the idea that the ion- 
acoustic waves may be produced by the ES decay of the L- waves (see also Cairns & Robinson 
1995). On the other hand, below 50 Hz fluctuations are most hkely of electromagnetic nature, 
such as whistlers. Summarizing, the analysis by Thejappa & MacDowall (1998) strongly 
suggests the presence of ion-acoustic waves in the range 50 — > 450 Hz, for this particular 
event, and that these waves are produced as the result of the ES decay of the beam-excited 
Langmuir waves. 

For the event under consideration, the numerical values of the relevant parameters are 
(Thejappa & MacDowall 1998): Vs = 5.1 x lO^cm/seg, Vtc = 1-6 x lO^cm/seg (corresponding 
to Tp = 1.7 X lO^K), Vsvi = 3.4 x lO^cm/seg, V(f, ~ Vhaan ± 25%, with a mean beam velocity 
estimated to be H)eam = 3.5 x lO^cm/seg, /pe = 2.3 x lO^'Hz, 6'Lr ~ 45°. 

Using these numerical values we compute now, in the two limit cases under analysis, the 
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predicted ion-acoustic wave frequencies. In all cases, we assume that the primary (beam- 
excited) Langmuir waves form a mean angle ^Lr = 45° with respect to the solar wind 
direction, and that the primary Langmuir waves have phase speeds in the range = ^eam± 



Figure 2 shows that, in the diffusive limit, decays proceed most likely for angles ^sl in the 
range 60° — > 75° (in this range the probability is about twice the back-scattering probability). 
Using this range of angles in equation (23), we find that the diffusive-scattering assumption 
yields ion-acoustic waves in the range 50 Hz — ^ 215Hz. On the other hand, if back-scattering 
is assumed, wc have ^^sl = 0°, and equation (23) yields ion-acoustic frequencies in the range 
215 Hz 415 Hz. The difference between the predictions in both limits is due to the term 
(cos (^sl) — ^^jT-j, that is maximized in the back-scattering case as cos (^^sl) ~ +1- 



If the the primary Langmuir waves form a variable angle 0° < Oj^,. < 90° with re- 
spect to the solar wind direction, the frequency ranges predicted in both limits overlap each 
other, within the observed range. In this case, the back-scattering limit yields predicted 
ion-acoustic frequencies over the whole observed range (as already pointed out by Thejappa 
& MacDowall 1998). We thus find that the predictions, for both asymptotic cases are consis- 
tent with the observations. Taking mean values of the different parameters, we find out that 
for the back-scattering assumption the predicted frequencies are consistent with the higher 
frequency portion of the observations (see also Cairns & Robinson 1995). On the other hand, 
the diffusive-scattering limit yields predicted frequencies that are consistent with the lower 
frequency portion of the observed electrostatic fluctuations. 

Further empirical support arises from a recent observational work by Thejappa et al. 
(2003), where they analyze URAP data for another type HI burst, in a similar fashion to 
the case already described. This case corresponds to in-situ observations at much larger 
heliocentric distances, specifically at 5.2 AU, where the local ambient plasma parameters 
are very different (/pc 2 kHz, Tg ~ 8.8 x lO'^K) to those of the previous case. The beam 
velocity (and hence the primary Langmuir waves phase velocity) is estimated to be ~ 
1.5 X lO^cm/seg ± 30%. In this other case, the WFA-URAP instrument also detected electric 
field fluctuations correlated in time with Langmuir bursts. The electric waves that can be 
safely assumed to be of electrostatic nature (i.e. with no simultaneous detection of magnetic 
field fluctuations above the background level) are in the frequency range 20 — > 200 Hz. 

We repeated our probability calculations for this case, and find out that, in the diffusive 
hmit, decays proceed most likely for angles ^sl in the range 55° — > 73° (in this range the 
probability is about twice the back-scattering probability). Using this range of angles, in the 
diffusive limit, the predicted ion-acoustic frequencies result in the range OHz 60 Hz. On 
the other hand, under the back-scattering assumption the predicted frequencies are in the 
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range 55 Hz —> 105 Hz. 

5. Conclusions 

From the theoretical point of view, an analysis of the ES decay rate indicates that 
the process is dominated by two limiting cases: diffusive (small angle) scattering and back- 
scattering. We find that the decay rates are comparable for both cases, being the diffusive 
decay rates systematically larger than those of the back-scattering limit. Given the results of 
our quantitative comparative analysis, we believe that both limiting cases arc acting at the 
same time. Furthermore, we speculate the net effect of diffusive decay and back-scattering 
acting simultaneously is then that of a diffusive isotropization of the bcam-gcncratcd Lang- 
muir spectrum. The back-scattering generates a backwards directed spectrum (respect to 
the beam propagation direction), but then this secondary spectrum can also diffuse by small- 
angle decay. Thus, if the ES decay is assumed to be present in solar flare and type III solar 
radio burst scenarios, we believe that its capability to isotropize Langmuir waves can not be 
neglected. 

From an observational point of view, we refer to analyses of a type 111 bursts by The- 
jappa & MacDowall (1998) and Thejappa et al. (2003). They registered low frequency elec- 
trostatic fluctuations which are strongly correlated in time with impulsive Langmuir wave 
bursts. We find that the observed frequencies of these fluctuations are consistent with the 
predicted frequencies of ion-acoustic waves produced by the ES decay of the primary (beam- 
excited) Langmuir waves. Under the back-scattering assumption the predicted frequencies 
are consistent with the higher frequency portion of the observations. On the other hand, 
the diffusive-scattering limit yields predicted frequencies that are consistent with the lower 
frequency portion of the observed electrostatic fluctuations. The relative burst intensities 
at the two frequency ranges could then serve as a proxy for the relative effectiveness of the 
occurrence of the ES decay in both limits, that we anticipate to be comparable from our 
theoretical analysis. 

From this analysis, we speculate that both the back-scattering and diffusive limits of 
the ES decay may be relevant in type HI bursts (and presumably also in solar flare radio 
events). Also, the beam-generated Langmuir turbulence may become isotropic as the ES 
decay proceeds mainly in the two limit cases of diffusive scattering and back-scattering. We 
postpone for a future work the investigation of the coupling between the ES decay and the 
quasi-linear beam-turbulence interaction. We believe that the simultaneous consideration 
of both effects can bring HXR and radio emission predictions into a better agreement with 
observations. This is due to the fact that its inclusion will imply a limitation of the effective- 
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ness of the quasi-linear relaxation, yielding to lower Langmuir turbulence levels, and hence 
lower radio emissivity. 
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Fig. 1. — Momentum conservation in the electrostatic decay, k^i = \s.l2 + kg. Angles a and 
P are given by equations (4) and (5). 

Fig. 2. — Dimensionless decay rate P(ks), as a function of {ks, ^sl)- 

Fig. 3. — Model for the angular distribution of the Langmuir spectrum, with 6c = 25° and 
— 2°. The vertical lines indicate the angles 9c (full line), and 9c ± (dashed lines). 

Fig. 4. — Diffusive decay: dimensionless factor Pi{ks,9s) for 9c — 5°, 25°, and A^ = 
2°. Profiles are ploted as a function of 9s, with different curves corresponding to kg — 
0, .., {ki) /lO (curves with larger Pi values correspond to larger ks). 

Fig. 5. — Diffusive decay: dimensionless factor Pi{ks,9s) for 9c = 0.5° and A9 = 0.5°. Pro- 
files are ploted as a function of 9s, with different curves corresponding to ks — 0, .., (/c^) /lO 
(curves with larger Pi values correspond to larger ks)- 

Fig. 6. — Diffusive decay: for a fixed ks = {ki) /lO, dimensionless factor Pi{ks,9s) for 
9c = 5°, 25°, and A9 = 2°, 10°. Profiles are ploted as a function of 9s, with thin curve 
corresponding to A^ = 2° and thick curve to A^ = 10°. 

Fig. 7. — Back-scattering: dimensionless factor Pi{ks,9s) for 9c = 5°, 25°, and A^^ = 2°. 
Profiles are ploted as a function of 9s, with thin curves corresponding to ks = 'ik^min — 2A;o 
and thick curves to ks = "^k^max ~ 2/co- 

Fig. 8. — Back-scattering: dimensionless factor Pi{ks,9s) for 9c — 0.5° and A^ = 0.5°. 
Profiles are ploted as a function of 9s, with thin curves corresponding to ks — 2kL^rnin ~ 2A;o 
and thick curves to ks — '^kL,max — 2A;o- 
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Rates for =5°, A0 = 2°, 10° 
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Rates for 0,, = : 
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